Resonant tunneling in a schematic model 
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Tunneling of an harmonically bound two-body system through an external Gaussian barrier is 
studied in a schematic model which allows for a better understanding of intricate quantum phe- 
nomena. The role of finite size and internal structure is investigated in a consistent treatment. The 
excitation of internal degrees of freedom gives rise to a peaked structure in the penetration fac- 
tor. The model results indicate that for soft systems the adiabatic limit is not necessarily reached 
although often assumed in fusion of nuclei and in electron screening effects at astrophysical energies. 
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The phenomenon of quantum tunneling is relevant in 
several areas of physics, from chemical reactions, elec- 
tronic circuits, to nuclear fission and fusion processes be- 
low the Coulomb barrier. In many cases one has to deal 
with the tunneling of composite objects. In general one 
tries to identify macroscopic degrees of freedom (denoted 
by R) and intrinsic ones (denoted by ^), decomposing the 
Hamiltonian accordingly: 



(1) 



The macroscopic part, Hmac{R), is the Hamiltonian for 
the (few) collective variables and contains besides the col- 
lective kinetic energy the macroscopic potential that has 
to be transversed by quantum tunneling. Hint{S.) con- 
tains the (many) intrinsic degrees of freedom of the many- 
body system, while V{R,£,) describes an interaction be- 
tween intrinsic and macroscopic variables. For example, 
if one wants to describe electron screening effects at astro- 
physical energies 1] , B could be the distance between the 
nuclei of two fusing atoms and ^ could denote the electron 
degrees of freedom. Or in nuclear physics R may denote 
the shape degrees of freedom for a heavy nucleus that un- 
dergoes spontaneous fission and ^ the individual nucleon 
degrees of freedom (see for example _2j). The question 
of how the tunneling reaction is influenced by finite size 
and structure of the composite object is an intriguing 
one. An important aspect is to determine which degree 
of freedom must be taken into account in a theoretical 
description and which can be neglected. The description 
of the tunneling of a composite object with many degrees 
of freedom is a very complex many-body problem. There- 
fore, very often radical approximations are made. As an 
example, nuclear cross sections at very low energy which 
are important in astrophysics and where no experimen- 
tal data exist, are often determined by extrapolations 
based on the one-dimensional result for point-like par- 
ticle tunneling. On the other hand, in heavy-ion fusion 
reactions it is well known that the coupling of the relative 
motion of the colliding nuclei with the nuclear intrinsic 
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motion strongly enhances the fusion cross section in the 
sub-barrier energy region. This was proven in a number 
of precise experiments, where the fusion cross section for 
intermediate mass systems (see e.g. 0) was found to be 
higher than simple one-dimensional prediction for tun- 
neling through a potential barrier formed by the attrac- 
tive nuclear interaction and the repulsive Coulomb force. 
The coupling between macroscopic and microscopic coor- 
dinates has been studied in multidimensional approaches 
(with many internal degrees of freedom), with different 
coupling schemes and with different approximations (see 
4] and references therein). 

In this paper we investigate the relation between the 
translational motion and the internal degrees of freedom 
of a composite object in a schematic but fully consistent 
model. In spite of the simplicity of the model it has all 
the ingredients to understand for example the sub-barrier 
fusion of soft nuclei, that can easily vibrate. In fact the 
low energy fusion cross section is usually dominated by 
s-wave fusion so that one deals only with one collective 
variable, the radial coordinate R, and a few low lying 
excited states |^ , like in the schematic model we describe 
in the following. 

The Hamiltonian of two interacting particles with iden- 
tical mass m under the influence of an external potential 
barrier is 
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where xi^2 denote the coordinate of particle 1 and 2, re- 
spectively, while Vint and V are the intrinsic potential 
and the external barrier, felt by both particles. Perform- 
ing a transformation to the center of mass (CM.) and 
relative coordinate denoted by R and ^, respectively, the 
Hamiltonian becomes 
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+ VintiO is the Hamiltonian of 



Here Hint{C) = 
the intrinsic system, while M = 2m and ^ = m/2 rep- 
resent the center of mass and the reduced mass, respec- 
tively. As one can see, the external potential, depending 
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on its functional form, may generate coupling terms of 
different orders between the macroscopic and the inter- 
nal coordinates. 

The standard theoretical approach to study the effect 
of internal excitations induced by the coupling potential 
is to solve the coupled-channel equations, which in our 
case read 
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where E is the total energy of the system and e„ is the 
internal excitation energy. These equations are obtained 
by introducing the eigenstates of the internal systems, 
i.e. 



and expanding the total wave function as 

N 
n=0 



(5) 
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N being the number of internal excitations considered. 
The expansion coefficients (pni{R) depend on the CM. 
coordinate. Here the sub-index i denotes the initial chan- 
nel, i.e. the starting energy level of the internal system. 
In Eq. (0J the potential matrix elements 
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can be interpreted as the effective potentials felt by the 
CM. due to presence of internal degrees of freedom. 

Equation 10} consists of a set of N coupled second or- 
der differential equations and, in case the particle is inci- 
dent on the barrier from the left hand side, the boundary 
conditions we require for its solution are 
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Here fc„ = ^ 2M(£' - e„)/?i is the wave number of the 
n-th channel. The inclusive penetration factor, or total 
transmission coefficient, is then given by 
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where the transmission probability for each channel is 
defined as 



ni — , \-^ni 
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In the schematic model that we would like to solve 
we assume an harmonic internal potential Vint{C) = 



^/xJl^^^ — iTifi and a Gaussian external barrier V{x) = 

2 2 

Vbe"" ^ . In this case the internal excitation energy be- 
comes £„ — nhVL and the effective potential matrix ele- 
ments have the following analytical form 5, f| 
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if n + j is even, and V,„(i?) = otherwise, where 
P = {nH./hy^'^ is the harmonic oscillator parameter. We 
would like emphasize that in this case the potential ma- 
trix elements can be calculated exactly and treated con- 
sistently within the model: they are given by the prod- 
uct of a Gaussian times a linear combination of Hermite 
polynomials. Therefore, they may show some structure 
depending on the chosen parameters. This is not equiva- 
lent to making a separable ansatz for the coupling poten- 
tial and then assuming constant form factors or a Gaus- 
sian parameterization as proposed in Ref. 7]. One can 
note that the potential matrix elements of connect- 
ing channel j with channel n depend on the internal sys- 
tem via the (3 parameter. The first matrix element, Vqo, 
which accounts for the elastic transition from the ground 
state ( j = 0) to itself (n = 0) , is a simple Gaussian func- 
tion whose width depends on the original width of the 
external barrier ('^ 1/^^) and on the harmonic oscillator 
parameter (3. In case of a very stiff internal harmonic os- 
cillator, i.e. for n ^ oo, one sees that Voo(-R) — > 2V{R). 
This is the limit case in which the composite object be- 
haves like a point particle, feeling two times the external 
potential. Thus, for finite CI the diagonal matrix ele- 
ments Vnn account for finite-size effect, while the cou- 
pling terms {n ^ j) are responsible for transition to the 
excited states. Performing a Taylor expansion of the cou- 
pling potential for ^ <C l/i^ up to the second order one 
gets 
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In case the two-body system is composed of identical par- 
ticles, both interacting with a symmetrical external bar- 
rier, the coupling potential does not include any linear 
coupling term as V'{R)£,. Actually, it consists in an ex- 
pansion over all even coupling terms. This means that 
the barrier effect on the intrinsic system is reduced, up to 
the second order of the expansion, to an additional har- 
monic potential that shrinks or stretches it, depending 
whether the second derivative of the external potential 
is positive or negative, respectively. By comparing the 
intrinsic harmonic potential Vint with the induced har- 

^2 

monic potential V"{R)^, for example at i? = where 
V"{R) presents a minimum with maximal amplitude, one 
can say that if hfl > fivyJVoJJl (or hVl < hi^^/Vo/fi) the 
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FIG. 1: Effective potential matrix elements for a two level 
system in comparison with the external barrier (parameters 
are: hn = 1 MeV, m = 938 MeV, Vo = 10 MeV and u = 0.1 
fm-i). 



intrinsic system is stiff (or soft) with respect to induced 
excitations. If the object is very stiff one expects the 
internal excitation to play a negligible role, whereas in 
case of a very soft system many internal levels have to be 
taken into account. 

In the following we solve the coupled channel problem 
using the potential matrix elements of Eq. Hll|l . In order 
to gain insight into the effect of the coupling we restrict 
ourselves to the case of a two level system. The internal 
degrees of freedom are taken to be initially in the ground 
state with excitation energy eo. Thus, in the presented 
model the only possible transition will be to the second 
excited state (n — 2), since the coupling term Voi{R) 
vanishes due to the mirror symmetry of the Gaussian 
barrier. For our two level system we can then define the 
potential matrix as following 
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where €2 = 2710 is the energy difference between the two 
levels, being eg = 0. Since we are interested in studying 
the effect of induced internal excitations, we will start 
to consider a soft object, where the internal structure 
plays a relevant role. In our first analysis we use follow- 
ing parameters: harmonic oscillator frequency hQ = 1 
MeV, mass m = 938 MeV, barrier amplitude Vq — 10 
MeV, barrier inverse width — 0.1 fm~^, which lead to 
hiyy/\^=2.88 MeV. 

In Fig. 1 (a) we firstly show the consistent potential 
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FIG. 2: Total transmission coefficient and separate T20 and 
Too contribution as a function of the total energy of the system 
for no, = 1 MeV. 



matrix elements of Eq. pi(l for the chosen parameters, 
in comparison with the external barrier. One can note 
that, due to finite size effect, the potentials Vqo and V22 
show a broadened structure with respect to the exter- 
nal potential, while the coupling term V02 changes sign 
two times. In Fig. 1 (b), we present the the so-called 
eigenbarriers (or eigenpotentials), obtained diagonalizing 
the potential matrix of Eq. p3|) at each position of the 
macroscopic coordinate R (see e.g. 0). 

We have then performed the coupled channel calcula- 
tion integrating Eq. ^ from R = —40 fm to i? = -40 
fm, and imposing the incoming wave boundary condition 
of Eq. ijHJi. It is known that with this method, often used 
in heavy-ion collisions calculations 0, it is sometimes 
difficult to obtain a stable solution with a controlled nu- 
merical accuracy. In order to check the numerical results 
we also solve the problem with the more stable method of 
the variable reflection amplitude , where the set of two 
coupled second order linear differential equations is trans- 
formed into a set of four coupled non linear first order 
differential equations. Other similar methods have been 
proposed to stabilize the numerical solutions of the cou- 
pled channel problem (see e.g. ^3)- We have obtained 
the same result with the above mentioned two methods 
with a relative percentage error of about 1-3 % in the 
presented energy region, indicating that the numerical 
accuracy is under control in the considered example. 

In Fig. 2 we show the result of the total transmission 
coefficient T and of the separate contributions T20 and 
Too. As one can see, the penetration factor presents 
a peaked structure at a total energy value of about 
^-pcak — 14 MeV. The pronounced peak is mainly due 
to excitation to the energy level £2, as can be seen from 
the fact that also T20 presents a peak, while the elas- 
tic channel Too is rather flat in this energy region. This 
agrees with the intuitive picture that the system avoids 
the higher barrier Vqo by going to the excited state which 
then tunnels easily through the lower barrier V22 + £2 (see 
Fig. 1). However, this picture does not hold for higher 
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FIG. 3: Transmission coefficients for different vaiues of the 
liarmonic osciliator frequency as a function of tlie total energy. 



energies where the situation is inverted due to the effect 
of the coupUng and finally, for energies somewhat higher 
than the external barrier (20 MeV), the elastic channel 
Too dominates the total transmission. Thus, a proper 
treatment of the coupling V02 is important. 

We would like to point out that the consistent treat- 
ment of the internal degrees of freedom in Eq. ((TJ pro- 
duces different widths for the diagonal potential matrix 
elements and oscillations in the coupling term. This leads 
to the emergence of a resonant structure in the penetra- 
tion factor. A peak in the transmission coefficient is not 
found if one parameterizes each potential matrix element 
with a simple Gaussian function with the same width 
as done in [3, H, > where a shoulder-like structure is 
found. 

We can then investigate the effect of the different stiff- 
ness of the composite object by keeping the same param- 
eters for the external barrier and changing the internal 
frequency of the harmonic oscillator. In Fig. 3 we show 
the transmission coefficients for different values of hft as 
a function of the total energy. One can see the the peaked 
structure moves towards higher energies with growin g in- 
trinsic frequency. When hft is larger than Tiv^/Vq/Jj. , the 
structure in T becomes less pronounced and more simi- 
lar to a shoulder, as found in j^. The reason is that, as 

increases, the widths of the diagonal matrix elements 
Vqo and V22 become more similar and the coupling V02 
gets smaller, leading to a potential matrix similar to the 
parameterization of Ref. 8] for high excitation energies. 
For = 2 MeV and higher the transmission coefficient 
at low energy is dominated by the elastic channel Too, 
since the consistent treatment of the potential matrix el- 
ements gives a Vbo barrier lower than the V22 + £2, in 
contradistinction to the case of Ml — 1 MeV depicted in 
Fig. 1. 

We also investigated the model with two excited in- 
trinsic states. In general the three level system picture 
within this consistent model leads to the emergence of 
two peaked structures: the position of the first peak is 
shifted towards lower energies and that of the second to- 



FIG. 4: Effects of the internal degrees of freedom for fiQ. = 1 
MeV: comparison of the coupled channel result with the no 
coupling calculation. Transmission coefficients for tunneling 
through the eigenbarriers are also shown. 



wards higher energies with respect to the location of the 
single peak found in the two level system. Such a be- 
havior was also found for the shoulder-like structure of 
the three level system in with respect the two level 
system in In case of three levels the structure of the 
consistent matrix elements becomes more complicated, 
therefore we prefer to stick to the simpler two level ap- 
proach to understand further the effect of the coupling. 

In order to estimate the dynamical effect due to ex- 
citation of internal degrees of freedom one should com- 
pare the coupled channel calculation with the uncoupled 
problem, where only the potential Vbo is considered. As 
already mentioned, this matrix element does not allow 
for internal transitions, but still accounts for finite size 
effects. This is different from the tunneling through the 
bare barrier 2V , since this would be equivalent to reduce 
the problem to the tunneling of a point-like particle. In 
Fig. 4 we show the above mentioned transmission coeffi- 
cients in case of a soft object, i.e. for Ml = 1 MeV. In the 
peak region the effect of the excitation of the internal de- 
grees of freedom leads to an enhancement of the transmis- 
sion coefficient T of about 5 orders of magnitude. On the 
other hand at higher energies the tunneling probability 
is lower than in case of no coupling. The enhancement of 
T in the coupled channel calculation with respect to the 
no coupling case strongly depends on the harmonic os- 
cillator frequency: it decreases with growing fi, i.e. with 
increasing stiffness of the composite object. In Fig. 4 we 
also compare our result with the transmission coefficients 
found for the tunneling through the eigenbarriers of the 
two-level system. In the following we recall the meaning 
of the eigenbarriers. 

Considering the part of the Hamiltonian of Eq. (|3Jl that 
depends on the intrinsic coordinate, i.e. 

H{R, = H,r.t{0 + v(r + IV v(r - Pi (14) 
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one can define the adiabatic states by minimizing it with 
respect to the internal degrees of freedom at each position 
R This translates then into an eigenvalue problem, 
which in our case reads 
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with W{R) as defined in Eq. ifT^ . The solution of the 
eigenvalue problem is found by diagonalizing the poten- 
tial matrix, i.e. by considering the eigenbarriers. In lit- 
erature the tunneling through the lowest eigenbarrier is 
often called adiabatic transition. 

As one can notice from Fig. 4, the adiabatic picture 
is very different from the result of the coupled channel 
calculation, though a small structure is found in the adi- 
abatic tunneling in the same position as the pronounced 
peak of the coupled calculation. Furthermore, we ob- 
serve that the adiabatic transmission coefficient is always 
larger than the coupled result. For higher energy the cou- 
pled channel calculation agrees with the result obtained 
for the tunneling through the second eigenbarrier. 

The concept of eigenbarriers is useful in case one wants 
to describe the fusion cross section as given by an average 
over the contribution form each eigenbarrier with appro- 
priate weight factors. A method to extract the barrier 
distribution from the measured cross section was pro- 
posed in Ref. From a purely theoretical point of 
view, the barrier distribution picture is correct only when 
the transformation that diagonalizes the matrix W{R) 
does not depend on the coordinate R. An approximation 
which is often made consists in considering the eigen- 
barriers and then evaluating the weight factor at a fixed 
position of R, assuming them to be constant. In Ref. 8] 
it was shown that in case one parameterizes the poten- 
tial matrix elements by Gaussians of the same width the 
weight factors are almost constant as a function of the 
energy. It was also proven that they are very different 
from those approximately estimated at the location of 
the maximum of the bare barrier For a two level 
system the weight factors are given by 

w^iE) = [TiE)-T,{E)]/[{T+{E)-T^iE)], 
w+{E) = [T+{E)~T{E)]/[{T+{E)-T^{E)], (16) 

where T-^^{E) denote the transmission coefficient for the 
first eigenbarrier and for the second, respectively, and 
T{E) is the total transmission as calculated in the cou- 
pled channel case. In Fig. 5 we show these optimum 
weight factors in our consistent model for the case of 
hil — 1 MeV. One can note that they are not constant 
in the energy region where the peaked structure in T is 
presented, showing that the approximation of constant 
weight factors does not hold. This is related to the fact 
that the transition is not adiabatic, as already discussed. 
We observe that with increasing stiffness of the composite 
object the weight factor show a less pronounced variation 
as a function of the energy. 

In order to gain more insight into the meaning of the 
adiabatic picture, we recall that in case of a two level 
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FIG. 5: Weight factors for the two eigenbarriers in case of 
hfl = 1 MeV. 



system the orthogonal rotation that is needed to diago- 
nalize the symmetric matrix W{R) can be parameterized 
by one single mixing angle 0{R) as follows 



cos 61 (i?) sin6'(i?) 
-sin6'(_R) cos 6'(i?) 



(17) 



At this point, if one would like to transform the coupled 
channel equation of , one should accordingly transform 
also the kinetic energy, which does not commute with the 
rotation operator, that depends on the coordinate R. In 
fact, denoting with K the diagonal kinetic energy matrix, 
the transformed matrix looks like 
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-im+iO'iR))' 
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(18) 

Thus, one can see that the adiabatic picture, which con- 
sists in considering only the first eigenbarrier solving an 
uncoupled problem, is valid in case that the first and 
second derivative of the mixing angle are negligible, i.e. 
in case that the transformed kinetic energy tends to the 
original diagonal one. In order to see whether this is the 
case or not in the considered example we show in Fig. 6 
the mixing angle (a) and its derivative (b) for different 
values of the internal frequency as a function of the 
CM. coordinate. Firstly one notes that in all cases 9 
crosses zero two times; these are the two points in which 
the coupling matrix element Vo2{R) = and thus W{R) 
is already diagonal. This is not the case if one param- 
eterizes all matrix elements with Gaussian functions of 
the same width as was done in 0, Q j where the mixing 
angle is always positive and can be maximal 45° in case 
of strong coupling. If one then focuses on the softer case 
of = 1 MeV, which corresponds to the example shown 
in Fig. 2 and Fig. 4, one observes that the angle changes 
suddenly sign, going from about 12° to about —64°. The 
derivative of the mixing angle is thus very large in the 
vicinity of the two zeros of 9. For this reason its contri- 
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FIG. 6: Mixing angle and its first derivative as a function of 
CM. coordinate for different values of h^l. 



bution in the transformed kinetic energy is not negligi- 
ble. In fact, the extra diagonal term ^{O'^axi^^ which 
is about 3.25 MeV at i? = ±6.6 fm, has to be added to 
the adiabatic potential, which is maximally 13.67 MeV. 
Moreover, the off-diagonal matrix elements contain the 
second derivative of 9{R) which gives maximum values 

of j^d'max ~ 4.34 MeV, invalidating the adiabatic as- 
sumption. 

The rapid change of the mixing angle is related to a 
so-called Landau-Zener pseudo-crossing of the two levels 
[T^. Denoting the rotated state with 



h{R) 



(19) 



and considering the optimal case in which the mixing 
angle varies suddenly as 6{Ri) — > 0(i?2) = d{Ri) ± 7''/2 
going from coordinate i?i to R2, one can note that 
the two rotated states invert each other in the sense 
^o(i?i) 02(^2) and (j>2iRi) (f)o{R2)- In other words, 
the first energy level suddenly becomes the second and 
vice versa. In case of a soft object the variation of the 
mixing angle does not give rise to the maximal pseudo- 
crossing, but still, for example, at i?i « — 14 fm we have 
e = 12° and MRi) = (O.98(/)o+O.2102) and at i?2 = we 
have 9 = -64° with 02(i?2) = (0.90(/)o -I- 0.44(/)2), so that 
4'2{Ri) and 4>2{R2) are very similar. Therefore, there is 
partial pseudo-crossing which is related to a non adia- 
batic transition, as discussed above. In fact, looking at 
Fig. 6 one can see that if the composite object becomes 
stiffer, then the mixing angle is smaller and does not vary 
so rapidly as a function of the CM. coordinate, i.e. there 
is no Landau-Zener crossing any more. For this reason 
one expects the adiabatic limit to be recovered for a stiff 
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FIG. 7; Comparison of the coupled channel result with the 
no coupling calculation and with the adiabatic transition for 
ftfi = 3 MeV. 



object. To illustrate that we show in Fig. 7 the transmis- 
sion coefficient obtained from the coupled channel cal- 
culation, the no coupling case and the adiabatic one for 
Ml = 3 MeV. One can note that for energies lower than 
the external barrier the adiabatic limit is now recovered 
(dotted curve lies on top of the solid line). This time, 

w(CaJ' « 0.018 MeV and i^Ca. « 0.19 MeV and 
thus the changes in the kinetic energy matrix H18|l are 
negligible. The coupling to the excited state still leads to 
a higher penetration factor with respect to the uncoupled 
case. For even stiffer systems, for example = 5 MeV, 
we observe that the coupled calculation finally coincides 
also with the uncoupled and the adiabatic one. We do not 
obtain an energy shift of the coupled result with respect 
to the uncoupled case as in |^, since in our consistent 
treatment the excitation energy and the potential ma- 
trix element are not independent from each other: with 
growing f2 the excitation energy increases but the cou- 
pling potential Vo2{R) goes to zero. When the object is 
too stiff to be excited no structure is found in the trans- 
mission coefficient: the internal system starts from the 
ground state and emerges still in the ground state at the 
end of the barrier. In this limiting case a coupled channel 
calculation is not necessary, since the internal degrees of 
freedom do not play any role. 

To summarize, we present a schematic model to de- 
scribe the tunnel effect of a two-body system in a two 
level approximation. Assuming an harmonic oscillator 
as internal interaction and a Gaussian external barrier 
one can give a consistent description of the potential felt 
by the macroscopic coordinate due to presence of inter- 
nal degrees of freedom. No a priori parameterization of 
the potential matrix elements is assumed. Therefore, in 
the coupled channel picture the dynamics of the internal 
degrees of freedom and their interaction with the exter- 
nal barrier are treated consistently. A stiffer system not 
only has a larger intrinsic excitation energy but also the 
potential matrix elements change accordingly. The cou- 
pled channel calculation shows a peaked structure in the 
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transmission coefficient, that accounts for the excitation 
of the intrinsic system. We show that for a soft object the 
mixing angle is large and changes rapidly sign, so that 
the adiabatic limit is not approached at low energies. The 
resonant transition to the excited state is explained by a 
Landau-Zener level- crossing in a non adiabatic picture. 



As expected, the adiabatic limit is recovered in case of 
a stiff object, where the energy is not sufficient to ex- 
cite the internal structure. The model results suggest to 
investigate more carefully fusion of soft nuclei or elec- 
tron screening at astrophysical energies, where often the 
adiabatic approximation is used. 
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